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Abstract 

Understanding the patterns and drivers of carbon isotope discrimination in C 4 species is critical for predicting the 
effects of global change on C 3 /C 4 ratio of plant community and consequently on ecosystem functioning and 
services. Cleistogenes squarrosa (Trin.) Keng is a dominant C 4 perennial bunchgrass of arid and semi-arid 
ecosystems across the Mongolian plateau of the Eurasian steppe. Its carbon isotope discrimination ( 3 A) during 
photosynthesis is relatively large among C 4 species and it is variable. Here the A of C. squarrosa and its potential 
drivers at a nested set of scales were examined. Within cohorts of tillers, 13 A of leaves increased from 5.1%„ to 8.1% 0 
from old to young leaves. At the local scale, "A of mature leaves varied from 5.8 , „ to 8.4' ,, increasing with 
decreasing grazing intensity. At the catchment scale, "a of mature leaves varied from 6.2 ,,„ to 8.5% 0 and increased 
with topsoil silt content. At the regional scale, "A of mature leaves varied from 5.5% 0 to 8.9 '„„ increasing with 
growing-season precipitation. At all scales, 1 A decreased with increasing leaf nitrogen content (N| ea f)- N lea f was 
positively correlated with grazing intensity and leaf position along tillers, but negatively correlated with precipitation. 
The presence of the correlations across a range of different environmental contexts strongly implicates N| eaf as 
a major driver of "a in C. squarrosa and, possibly, other C 4 species. 
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Introduction 

Grasslands have been experiencing rapid shifts in botanical 
composition, structure, and functioning driven primarily by 
global climate change and human disturbances (Alward 
et al, 1999; Grime et al, 2000; Sage and Kubien, 2003). 
Such changes may also shift the dominance of the two 
distinct photosynthetic pathways, C 3 and C 4 , which differ in 
the mode of initial C0 2 fixation and have profound effects 
on ecosystem processes (Farquhar et al, 1989; Hunt et al, 
1996; Sage and Kubien, 2003; Kuzyakov et al, 2006). Most 
of the world's C 4 species are grasses (Teeri and Stowe, 1976; 



Ehleringer et al, 1997). The ratio of C 3 to C 4 photosyn- 
thetic types is an important component in understanding 
ecosystem processes because it can affect biomass pro- 
duction, carbon sequestration, and water cycling (Tieszen 
et al, 1997; Pyankov et al, 2000; Sage and Kubien, 2003). 
It can also influence the terrestrial carbon isotope signal, 
which is used for the land-ocean partitioning of C0 2 fluxes 
(Still et al, 2003; Suits et al, 2005). 

The actual C 3 /C 4 ratio is often obtained from the stable 
carbon isotope compositions of C 3 and C 4 plants and bulk 



Abbreviations: 13 A, carbon isotope discrimination; MAP, mean annual precipitation; MAT, mean annual temperature; MJulT, mean July temperature; Nieaf, leaf nitrogen 
concentration; PEPc, phosphoenolpyruvate carboxylase. 
© 201 1 The Author(s). 

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by- 
nc/2.5), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. 



4144 | Yang et al. 



samples in a two-member mixing model (Murphy and 
Bowman, 2007), where bulk samples may comprise soil 
organic matter (Wittmer et al, 2010), «-alkanes (Zhang 
et al, 2003), wool (Auerswald et al, 2009), or other 
materials. This approach is based on the difference in 
discrimination against the heavy carbon stable isotope ( A) 
during photosynthesis between C 3 and C 4 plants (Farquhar 
et al, 1989). A single mean C 4 end-member value is usually 
assumed in the two-member mixing model (Murphy and 
Bowman, 2007; Wittmer et al, 2010), whereas variable C 3 
end-member values are estimated, based on the well known 
correlations between the A of C 3 plants and environmental 
conditions (Mole et al, 1994; Mannel et al, 2007; Wittmer 
et al, 2008). The reason for this is that A of C 4 plants is 
reported to be less responsive to environmental variables 
than that of C 3 plants (Van de Water et al , 2002; Murphy 
and Bowman, 2009). However, this view has been challenged 
by a number of studies that have examined variations in A 
of C 4 plants, and which suggest that the A of C 4 species can 
vary with environmental variables (Henderson et al, 1992; 
Buchmann et al, 1996; Wang et al, 2005) and that this 
variability is related to ecophysiological responses and re- 
source efficiencies (Henderson et al, 1992; Ghannoum 
et al, 2002; Wang et al, 2005). For instance, the A of 
Amaranthus emeritus was ~6% c higher in low-light leaves 
than in high-light leaves (Tazoe et al, 2006). The A of 
Bothriochloa ischaemum increased 2%„ with an increase in 
annual precipitation from 350 mm year -1 to 700 mm year - 
(Liu et al, 2005). In Saccharum spp., A declined by 2% Q 
with increasing leaf nitrogen content (Ni eaf ) in a nitrogen- 
stress experiment (Meinzer and Zhu, 1998). 

The A of a plant can be obtained from the carbon 
isotope composition of a plant (§ 13 C p i ant ) and that of the air 
(§ 13 C air ) in which it grows: 

cl3 r cl3 r 
13 . 0 ^-air 0 ^-plant 



13 A of C 4 plant biomass is influenced by many factors, 
including isotope effects during diffusion of CO2 through 
the boundary layer, stomatal pore, cell walls, membranes, 
and in the liquid phase, dissolution, hydration by carbonic 
anhydrase, fixation of bicarbonate by phosphenolpyruvate 
carboxylase (PEPc) in mesophyll cells, fixation of C0 2 by 
Rubisco in bundle sheath cells, leakage of C0 2 from bundle 
sheath cells to mesophyll cells, and carbon allocation and 
partitioning, including photorespiration and dark respira- 
tion (Farquhar, 1983; von Caemmerer and Furbank, 2003). 
This complexity has not been analysed in full detail, mainly 
due to lack of instrumentation and experimental protocols 
for the measurements of all the necessary variables. To 
solve the problem, Farquhar (1983) proposed a simplified, 
but testable version of the model, which was subsequently 
used widely to investigate the influence of the supposedly 
strongest drivers of variation in A. This model states that 
A of C 4 plants depends on bundle sheath leakiness (<&, the 
proportion of C fixed by PEP carboxylation, which sub- 



sequently leaks out of the bundle sheath) and the ratio of 
internal and ambient C0 2 concentrations (C ; /C a ): 

13 A = a + (b 4 + h x <P - a) x Q I C a (2) 

where a is the discrimination of C during diffusion of 
C0 2 through stomata (4.4% 0 ), b 3 is the discrimination during 
fixation by Rubisco [27% 0 for C 4 plants; von Caemmerer 
et al (2008)], and b 4 is the discrimination during hydration 
of C0 2 to HCO 3 and fixation by PEPc. 

Eqn (2), predicts that A variation during photosynthesis of 
C 4 plants depends on b 4 , 3>, and on C ; /C a . These variables in 
turn depend on a host of environmental conditions. Temper- 
ature is the biggest influence on b 4 , inducing variation from 
-4.8% 0 at 33.8 °C to -6.2% 0 at 21.0 °C (Henderson et al, 1992). 
$ is an important ecophysiological trait that is directly 
associated with photosynthetic quantum yield and nitrogen 
use efficiency of C 4 plants (Henderson et al, 1992; Ghannoum 
et al, 2002; Wang et al, 2005). Any C0 2 that leaks from 
bundle sheath cells decreases the energetic efficiency of C 4 
photosynthesis, through futile costs of ATP for the regenera- 
tion of PEP (Furbank et al, 1990). <5 depends on the balance 
of the biochemical capacities of the C 4 cycle in the 
mesophyll cells and the C 3 cycle in bundle sheath cells, and 
the C0 2 conductance of the bundle sheath exterior wall (von 
Caemmerer and Furbank, 2003). Factors that could affect <& 
include the activity ratio of Rubisco to PEPc (Farquhar, 1983), 
such as N leaf (Sage et al, 1987), temperature (Pittermann and 
Sage, 2001), and light (Tazoe et al, 2008). Q/C a depends 
on the ratio of stomatal conductance to photosynthetic activ- 
ity and, hence, influences intrinsic water use efficiency 
(Ghannoum et al, 2002). Because there is an interaction 
between the influences of 3> and Q/C a on A, A can increase 
or decrease with Cj/C a depending on the term (b 4 +b 3 <&-a) 
(Farquhar, 1983). For 3»0.37, this term is greater than zero 
with 13 A>4.4% 0 and it increases with Q/C a . For O<0.37, this 
term is smaller than zero with "A<4.4% 0 and it decreases with 
increasing Q/C a . For A differing from 4.4% 0 , A will thus 
vary with C;/C a if $ is constant. 

Studies of the A of C 4 plants and its response to 
environmental variables have been conducted primarily in 
grasslands in North America, Europe, and Australia (Schulze 
et al, 1996; Murphy and Bowman, 2009; Ghannoum et al, 
2011) or in experiments under controlled conditions (Ranjith 
et al, 1995; Tazoe et al, 2006). Few studies of this type have 
been reported from the mixed C 3 /C 4 grassland of the 
Mongolian plateau, where the C 4 community of the semi- 
arid steppe is dominated by Cleistogenes squarrosa (Trin.) 
Keng, a C 4 (NAD-ME) perennial bunchgrass. C. squarrosa 
has a wide range occurring at meadow steppe, typical steppe, 
desert steppe, and sand dune ecosystems across the Mongo- 
lian plateau. The A in C. squarrosa is well above 4.4% c 
(Pyankov et al, 2000; Chen et al, 2002; Gong et al, 2008) 
and it is known to vary with soil moisture (Chen et al, 2002). 

In this study the A variability in C. squarrosa and its 
controlling factors were investigated across a range of 
ecological scales on the Mongolian plateau. It was hypothe- 
sized that A variability in C. squarrosa is driven by factors 
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controlling N leaf , which influence Q/C a and Factors 
influencing Ni eaf may include leaf age, grazing pressure, soil 
properties, or weather conditions, all of which vary at different 
scales. Hence the A of C. squarrosa was studied with a nested 
set of scales (i.e. tiller, local, catchment, and regional) to 
identify and quantify controls of Ni eaf and, in turn, A. 



Material and methods 

Study area 

The study area was located between 106.15 °E and 1 17.13 °E 
longitude and 40.62 °N and 45.57 °N latitude in the Inner 
Mongolia Autonomous Region, China (Fig. 1). This area of 
—200 000 km 2 constitutes the regional scale. Samples were 
taken at altitudes ranging from 1010 to 1635 m elevation. 
The temperate, semi-arid climate is characterized by warm 
summers and cold, dry winters (Chen, 1988) with mean 
annual temperature (MAT) between 0.3 °C and 5.6 °C and 
mean July temperature (MJulT) between 18.7 °C and 
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Fig. 1. Nested hierarchy of sampling sites in the Inner Mongolia 
grassland. Leaf position along a tiller constitutes the tiller scale, 
plots in a grazing experiment constitute the local scale, the Xilin 
river basin constitutes the catchment scale, and the steppe area of 
Inner Mongolia constitutes the regional scale. Numbers on the tiller 
scale are phytomer units from base. 



22.7 °C. Mean annual precipitation (MAP) increases from 
160 mm year -1 in the western part to 350 mm year -1 in the 
eastern part, and 60-80% of the annual total occurs in the 
May-August growing season (Bai et al, 2007). Regional 
sampling took place in five different years (2005-09) to take 
account of the large inter-annual variation in precipitation 
in the region. 

The Xilin river basin covers an area of 10 786 km 2 in the 
northeast of the region. Within this catchment, MAP 
increases from 275 to 346 mm year -1 along a northwest- 
southeast precipitation gradient, with elevation ranging 
from 1019 to 1482 m. MAT varies from 0.3 °C in the 
southeast to 2.5 °C in the northwest of the basin, with 
MJulT ranging from 18.7 °C to 20.9 °C. The soils in the 
Xilin river basin include phaeozems, shallow and stony soils 
at rocky outcrops, and sandy soils of low water-holding 
capacity in sand dune areas (Wiesmeier et al, 2011). 

A grazing experiment (116.67 °E, 43.56 °N) near the 
Inner Mongolia Grassland Ecosystem Research Station 
(IMGERS) of the Chinese Academy of Sciences at the 
middle reach of the Xilin river represents the local scale. 
MAT, MJulT, and MAP at IMGERS are 0.4 °C, 19.0 °C, 
and 337 mm year -1 , respectively (Zheng et al, 2011). 
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Growing season precipitation was geostatistically deter- 
mined according to Wittmer et al. (2008) for each sampling 
location and time of sampling. This was done by interpolat- 
ing between meteorological stations and taking into account 
1.5x1.5 km 2 resolution maps of MAP, mean monthly 
precipitation, MJulT, and MAT (Climate Source, Inc., 
Corvallis, OR, USA). The mean daily precipitation 
(mm d -1 ) between 1 May, which is the date that usually 
corresponds with the onset of growth of C. squarrosa (Liang 
et al, 2002), and the day of sampling was calculated. This 
measure of precipitation takes account of the fact that 
samples were collected on different days and in different 
years. 

§ 13 C a i r was obtained from a long-term monitoring station 
in Ulan Uul (Climate data online; NOAA NNDC, 2009), 
Mongolia, —460 km northwest of IMGERS. 



Soil analyses 

On the catchment scale sampling locations were selected to 
cover the full range of soils as known from soil mapping 
campaigns. Soil pits were dug to 1 m at locations represent- 
ing the mapping units and the soils were sampled horizon- 
wise. A comprehensive set of soil properties was measured 
(for details see Wiesmeier et al, 2011), including bulk 
density, texture, pH, and nitrogen content. Either topsoil 
properties (e.g. silt content) because water uptake by 
C. squarrosa is restricted to the topsoil (Yang et al, 2011), 
or whole-soil properties (e.g. N stocks) were used. 
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Grazing experiment 

The grazing experiment had 28 annually grazed plots 
(Fig. 1 , grazing experiment) within two experimental blocks 
differing in topography (flat and gently sloping). The plots 
had been grazed between June and September each year 
since 2005 at seven stocking rates: 0, 37.5, 75.0, 112.5, 
150.0, 187.5, and 225.0 sheep km -2 with three times higher 
stocking densities during the 4 months of grazing. Further 
information about the grazing experiment can be found 
elsewhere (Schonbach et al, 2009). 



Statistical analysis 

Linear regressions were used to evaluate relationships 
between A, Ni eaf , and the environmental variables of 
interest at a certain scale: leaf position, grazing intensity, 
soil properties, and growing season precipitation. The 
coefficient of determination was tested with a two-sided 
test for significance of the regression. Confidence inter- 
vals (95%) were calculated for the relationship between 
A and Ni eaf and for the population. All procedures were 
carried out in SPSS Version 16.0 (SPSS, Inc., Chicago, 
IL, USA). 



C. squarrosa sampling 

The tiller scale was sampled in an ungrazed plot, where no 
leaves were lost by grazing. Leaves from several individuals 
were sampled and pooled according to their phytomer 
position (Fig. 1). Leaves were collected in mid- July, end of 
July, and mid-August 2008. Resampling over time allowed 
us to distinguish between new leaves and old leaves grown 
between sampling dates. Old leaves at low phytomer 
positions were then compared between consecutive sam- 
pling dates to identify changes over time. 

At the local, catchment, and regional scales, mature, 
fully expanded, and sun-exposed leaves from several non- 
flowering individuals were sampled. At the local scale, leaves 
were collected from the plots of the grazing experiment in 
June, July, and August of 2007, and in May, June, and July 
of 2008. These two years differed in precipitation: 2007 was 
dry with 240 mm year -1 , and 2008 was an average year with 
362 mm year -1 . At the catchment and regional scales, 
samples were collected in each of the five years 2005-09. 
The precise locations and elevations of all sampling 
positions were recorded with a mobile GPS. In total, 51, 
144, 64, and 47 samples at the tiller, local, catchment, and 
regional scales, respectively, were collected. 



Isotope and N| eaf analysis 

All samples were dried at 60 °C for 48 h in a forced-draught 
oven and ground to homogeneity with a ball mill. Carbon 
isotope composition and nitrogen content were then mea- 
sured using an elemental analyser (NA 1110; Carlo Erba, 
Milan, Italy) interfaced (ConFlo III; Finnigan MAT, 
Bremen, Germany) to an isotope ratio mass spectrometer 
(Delta Plus; Finnigan MAT) following the protocol by 
Wittmer et al. (2008). Carbon isotope data were specified as 
8 C relative to the Vienna Pee Dee Belemnite standard: 

(5 13 C = Sample _ j ^ 
^sample 

13 1 T 

where R samp i e and ^ stan dard are the ratios of C/ C in the 
sample and standard. 

The precision for sample repeats was better than 0.1 5% G 
for 8 13 C and 0.04% for nitrogen content in dry matter. 



Results 

Relationship between 13 A and leaf position at the tiller 
scale 

The A of C. squarrosa decreased and Ni eaf increased from 
the lowermost leaf to the uppermost leaf (Fig. 2). Thus, 
young leaves had a low A but high N leaf . The A mainly 
changed from leaf 5 to younger leaves, while Ni ea f mainly 
changed among leaves older than leaf 7. Hence in a bulk 
sample of undefoliated plants, the variation in A can 
mainly be attributed to younger leaves whereas the varia- 
tion in Ni ea f can mainly be attributed to older leaves. 

The N leaf of a leaf produced at a given time remained 
unchanged between mid-July (average=1.6%, SD=0.2%) 
and August (average=1.6%, SD=0.3%), which suggested 
no nitrogen redistribution over time. N] eaf of the three 
uppermost (newly grown) leaves (average=2.4%, SD=0.1%) 
was significantly higher (7-"<0.01) in August than those 
grown in July (average=1.7%, SD=0.2%). 

13 A, stocking rate, and topography at local scale 

At the local scale/'A varied between 5.8% 0 and 8.4% c with 
an SD of 0.5^ o . N w varied between 1.8% and 3.8% with an 
SD of 0.4%. ' A decreased (i? 2 =0.04, #=144, i 5 <0.05) while 
the Nieaf increased (i? 2 =0.11, #=144, J P«0.001) with 
stocking rate (Fig. 3). A and N leaf also varied with year 
( J P<0.001) and month (P<0.001), but not with topographic 
position (P>0.05). 

13 A and soil variables at catchment scale 

Across the Xilin river catchment, *A of C. squarrosa varied 
between 6.2%„ and 8.5% G (SD=0.5% o ), N leaf varied between 
1.5% and 3.7% (SD=0.5%), and topsoil silt content varied 
from 2% to 40% (SD=10%). "A increased with topsoil silt 
content CR 2 =0.42, #=24, P<0.001) despite there being 
differences in slope between the two years (P<0.05), but 
Ni ea f was not correlated with silt content (P>0.05) (Fig. 4). 
Similar patterns were found with soil nitrogen stocks. The 
A increased significantly with soil nitrogen stocks 
(i? 2 =0.40, #=25, P<0M), but there was no relationship 
between soil nitrogen stocks and Ni ea f (P>0.05). The 
influence of soil nitrogen stocks and silt content on 
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Phytomer No. 

Fig. 2. Influence of leaf position (phytomer 1 is the first leaf from 
base) on A and N teaf of C. squarrosa. Means of three sampling 
dates in an ungrazed plot; error bars denote the standard error. 
Note that more than seven phytomers were found in August and 
only one sample was available for phytomer 9. 
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Fig. 3. Influence of stocking rate on A and N| eaf of C. squarrosa. 
Each data point shows the mean and standard error of 24 
samples taken from four plots with different topographic positions 
and management systems. 

"A could not be separated due to their correlation with each 
other; this also applied to other soil properties. 

13 A and precipitation at regional scale 

As expected, the range of "A of C. squarrosa (3.4% c ) was 
considerably larger at the regional scale than at the local 
(2.6% 0 ) and catchment scales (2.3% 0 ). Precipitation during 
the growing season had a significant, positive effect on 
"A (i? 2 =0.31, N=47, P<0.001) and a negative effect on N leaf 
OR 2 =0.56, N=47, P«0.001) (Fig. 5). The °A of C. 
squarrosa increased OP<0.05) whereas N leaf decreased 
(P<0.01) from the desert steppe to the typical steppe, and 
to the meadow steppe. 
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Fig. 4. Influence of topsoil silt content on 13 A and N| eaf of 

C. squarrosa at the catchment scale in 2007 (filled symbols), a dry 

year, and in 2008 (open symbols), a normal year. 
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Fig. 5. Influence of growing season precipitation on "A and N| eaf of 
C. squarrosa at the regional scale. Growing season precipitation was 
calculated based on year and time of sampling (see Materials and 
methods). 



13 A and N, eaf 

At tiller, local, and regional scales, 3 A was negatively 
correlated with N leaf (Fig. 6) although the range of values 
and factors that affected the variation in N leaf differed 
between scales (Figs 2, 3, 5). In spite of no significant 
relationship at the catchment scale (P>0.05), all except 
four of the total of 47 data points were included in the 95% 
confidence intervals for all pooled samples. Ni eaf explained 
23% of the variation in A when all samples were pooled. 
When regressed separately, the slope of the regression was 
significantly steeper at the tiller scale than other scales 
(P<0.01), but remained within the 95% confidence interval 
of the common regression slope for the other scales 
(Fig. 6). 
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Fig. 6. Correlations between A and N| eaf of mature leaves of 
C. squarrosa at the regional (steppe area of Inner Mongolia), 
catchment (Xilin river catchment), local (grazing experiment), and tiller 
scales (an ungrazed plot). The central line in each panel denotes the 
regression for the pooled samples, the inner envelopes (solid lines) 
are the 95% confidence intervals for the regression, and the outer 
envelopes (dashed lines) are the 95% confidence intervals for 
individual samples. 

Discussion 

This study represents the first comprehensive examination 
of the A variability of a C 4 (NAD-ME) perennial 
bunchgrass, C. squarrosa, and its controlling factors across 
a range of scales (i.e. tiller, local, catchment, and regional 
scale) on the Mongolian plateau. The results reported here 
have shown that A of C. squarrosa was large (on average: 
7.2% c ) and, at the region scale, it varied over a 4% c range 
(5-9% 0 ) across the Inner Mongolian steppe. This range is 
larger than that for other C 4 species in the study area. For 
instance, Salsola collina, a dicot sampled together with C. 
squarrosa in the grazing experiment, varied by 3% 0 under 
identical conditions (data not shown). A previous study has 
reported that B. ischaemum, a dominant C 4 perennial 
bunchgrass on the Chinese loess plateau, varied by only 
2% D across its range (Liu et al., 2005). This suggested that 
C. squarrosa could be a model species for studying the effect 



of environmental conditions on the A of C 4 plants, and 
that A of the C 4 community is likely to vary, largely due to 
the change in relative abundance of C. squarrosa. 

Many studies have demonstrated that foliar § 13 C of 
plants is negatively correlated with water availability (e.g. 
soil moisture and annual precipitation) across an environ- 
mental gradient (Farquhar et al, 1989; Stewart et al, 1995; 
Liu et al, 2005). The mechanisms underpinning these 
patterns, however, have yet to be explored. The results 
reported here indicate that N leaf has a strong and consistent 
effect on A of C. squarrosa at all scales and explained 23% 
of the variation in A, even though the factors influencing 
Ni ea f varied at different scales. Therefore confidence can be 
placed in this relationship, and it may be expected that 
other environmental variables that affect Ni eaf would also 
affect A in a similar way. 



Factors influencing N/ eaf 

Factors that affect N leaf include factors that control 
N availability, such as fertilization (Bai et al, 2010), factors 
that influence N uptake, such as root density and root 
activity (Farooq et al, 2009), and factors that relate to 
ontogenesis and organogenesis, such as the availability of 
water and nutrients, which determine the velocity of leaf ex- 
pansion and leaf dimensions where grazing could contribute 
(Anten et al, 1998). 

The influence of leaf position on Ni ea f may have been 
influenced by soil temperature, which was —12.6 °C in late 
May and increased to 21.1 °C in mid- August (average of 
10 d at 10 cm depth, data from IMGERS). N uptake of C 4 
is especially sensitive to low temperature (Farooq et al, 
2009). Remarkably, neither N] ea f nor A of newly formed 
leaves changed between mid-July and August, indicating 
that the influence of leaf position was not caused by 
redistribution of N from old to young leaves. The increase 
in leaf nitrogen with grazing intensity is well known and 
relates to the faster cycling of N and decrease in mean leaf 
age (Mattson, 1980). Decreasing N leaf with increasing 
precipitation is likely to be related to the dilution effect 
caused by additional plant growth. Temperature-related 
N uptake may also add to the effect of precipitation, as 
regional temperatures and precipitation are closely nega- 
tively correlated (R 2 —0.93, Bai et al, 2007) and hence their 
effects cannot be separated. Finally, the Birch effect of 
enhanced mineralization after a pronounced dry spell 
(Birch, 1964; Austin et al, 2004) could also contribute to 
improved N availability in drier areas where rainfall occurs 
infrequently. 

At the catchment scale, the influences on N leaf are less 
clear. A increased with silt content whereas Ni eaf did not, 
although the relationship between A and N leaf was similar 
for this scale as for the other scales. At this scale the 
variation of soil is expected to have a strong influence but, 
in contrast to leaf position, grazing intensity, or precipita- 
tion acting on the other scales, the soil cannot be arranged 
along only one single variable like silt content. Soil is an 
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aggregation of several soil properties, which may exert 
influence on N] ea f and in turn on A. 

The effect of A/ /eaf on 13 A 

In principle, the variability in 13 A could be derived from: 
(i) discrimination associated with photosynthesis; (ii) post- 
photosynthetic fractionation, which includes C0 2 release in 
dark respiration, or (iii) any other simplifications included 
in Eqn (2). Although the latter processes have not been 
studied in detail in C 4 plants (but see Kromdijk et al, 2008; 
Pengelly et cti, 2010), evidence for the contribution of post- 
nhotosynthetic fractionation comes from the fact that 
A obtained by online measurement of photosynthetic 
CC>2/ 12 C02 exchange is often less than A estimated from 
the carbon isotope composition of leaf dry mass (Evans 
et al, 1986; Kubasek et al, 2007). For instance, "A 
obtained from online measurement was — 0.75— 1 .5% 0 lower 
than that from leaf dry mass (Kubasek et al, 2007). 
Respiration and carbon allocation are responsive to envi- 
ronmental conditions including temperature and water 
stress (Ryan, 1991; Flexas et al, 2005) and thus can 
contribute to the observed variability. However, changing 
chemical composition seems unlikely to explain the correla- 
tion of A with Ni ea f. Proteins differ isotopically from 
cellulose only by 4% c (Benner et al, 1987). A difference in 
Ni ea f of 2%, corresponding to a difference in protein content 
of 12%, is insufficient to explain a difference in A of 3% 0 . 
Also, secondary compounds, which increase with leaf age 
and may thus relate to Ni ea f, cannot explain the variability 

Table 1. Variation in the ratio of internal to ambient C0 2 
concentration (C/Ca) and leakiness (3>) for the local scale and 
correlations with N| eaf 

Gas exchange measurements were taken from Zheng ef al. (201 1). 



Variable 


Minimum-maximum 


Correlation with N| eaf 


C/C a 


0.41-0.66 


-0.866 


<t> 


0.56-0.58 


+0.832 



of "A at the tiller scale, because the young leaves present in 
mid-July did not change in isotonic composition until mid- 
August. Hence variation in leaf A was most likely caused 
by the variation in discrimination during photosynthesis. 

In Eqn (2), A mainly depends on <I> and Ci/C a as 64 
changes little within the 4 °C range in air temperature in the 
study area. Ci/C a should increase with stomatal openness 
(Turner et al, 1984), but decrease with the activity of 
carboxylating enzymes (Sage et al, 1987; Sage and Pearcy, 
1987), which should increase with Ni ea f (Ranjith et al, 1995; 
Sage et al, 1987). $ depends on the activity ratio of 
Rubisco and PEPc, which is influenced by the allotment of 
nitrogen to them (Sage et al, 1987). At the local scale, Ci/C a 
was measured for young mature leaves close to the top of 
tillers by Zheng et al. (2011) simultaneously with our 
sampling. Ci/C a varied between 0.41 and 0.66 (Table 1), 
similar to the range of 0.34-0.59 found by Chen et al. (2005) 
(Table 2). Using these C ; /C a data and the A data reported 
here, <5 was estimated using Eqn (2). This indicated that $ 
ranged between 0.56 and 0.58 (Table 1). N] ea f from this 
study correlated with both C ; /C a and $ although variation 
of 0 was small. Therefore, A seemed to be influenced 
mostly by C;/C a , which is related to water use efficiency. 



Conclusions 

The findings reported herein demonstrate that N lea f is a key 
factor driving leaf A variation of C. squarrosa. Increasing 
Neaf was associated with decreased A, which seemed to be 
related to improved water use efficiency. The leaf-level 
response of A to N] ea f was found at all scales from tiller to 
local, catchment, and regional although the reasons for 
Neaf variation differed. Within cohorts of tillers, A de- 
creased with leaf position. Within the field experiment site, 
A decreased with grazing intensity. Within the Xilin river 
catchment, A decreased with decreasing topsoil silt con- 
tent. Across the regional study area, A decreased with 
decreasing growing season precipitation. Accordingly, it is 
proposed that water use efficiency increased for leaves 



Table 2. Leakiness (<&) and ratio of internal to ambient C0 2 concentrations (cyc a ) of NAD-ME plants from published sources 



Species 


O 


Ci/C a 


References 


Amaranthus caudatus 


0.37 


0.43-0.46 


Kubasek ef al. (2007) 


Amaranthus edulis 


0.32 


0.43 


Henderson ef al. (1 992) 


Amaranthus retroflexus 


0.40 


0.40-0.90 


Sage and Pearcy (1 987) 


Atriplex lentiformis 


0.54-0.62 




Zhuefa/. (1999) 


Atriplex rosea 


0.45 


0.47 


Henderson et al. (1 992) 


C. squarrosa 




0.34-0.59 


Chen ef al. (2005) 


Cynodon dactylon 




0.12-0.28 


Mantlana ef al. (2008) 


Eleusine coracana 


0.40 


0.38 


Henderson ef al. (1 992) 


Eleusine indica 


0.37, 0.27 


0.48, 0.58 


Farquhar (1983) 


Eragrostis lehmanniana 


0.45-0.57 


0.08-0.70 


Fravolini ef al. (2002), Mantlana ef al. (2008) 


Panicum schinzii 


0.32 


0.38 


Henderson et al. (1992) 


Total range 


0.27-0.62 


0.08-0.90 




Mean 


0.42 


0.44 
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produced later in the growing season, and that it also 
increased with grazing intensity. But it declined with 
increasing silt content and water-holding capacity of the 
soil and precipitation during the growing period. 
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